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Abstract. In this paper we present a new theory of calculus over 
/c-dimensional domains in a smooth n-manifold, unifying the dis- 
crete, exterior, and continuum theories. The calculus begins at a 
single point and is extended to chains of finitely many points by 
linearity, or superposition. It converges to the smooth continuum 
with respect to a norm on the space of "monopolar chains," culmi- 
nating in the chainlet complex. Through this complex, we discover 
a broad theory of coordinate free, multivector analysis in smooth 
manifolds for which both the classical Newtonian calculus and the 
Cartan exterior calculus become special cases. The chainlet opera- 
tors, products and integrals apply to both symmetric and antisym- 
metric tensor cochains. As corollaries, we obtain the full calculus 
on Euclidean space, cell complexes, bilayer structures (e.g., soap 
films) and nonsmooth domains, with equal ease. The power comes 
from the recently discovered prederivative and preintegral that are 
antecedent to the Newtonian theory. These lead to new models for 
the continuum of space and time, and permit analysis of domains 
that may not be locally Euclidean, or locally connected, or with 
locally finite mass. 



Preface 

We put forward a novel meaning of the real continuum which is found 
by first developing a full theory of calculus at a single point - the origin, 
say, of a vector space - then carrying it over to domains supported in 
finitely many points in an affine space, and finally extending it to the 
class of "chainlets" found by taking limits of the discrete theory with 
respect to a norm. Local Euclidean structure is not necessary for the 
calculus to hold. The calculus extends to /c-dimensional domains in 
n-manifolds. We do not rely on any results or definitions of classical 
calculus to develop our theory. In the appendix we show how to derive 
the standard results of single and multivariable calculus in Euclidean 
space as direct corollaries. 



Date: December 21,2005, edited May 12, 2006. 

1 



2 J. HARRISON DEPARTMENT OF MATHEMATICS U.C. BERKELEY 

This preprint is in draft form, sometimes rough. There are some 
details which are still linked to earlier versions of the theory. It is 
being expanded into a text which includes new applications, numerous 
examples, figures, exercises, and necessary background beyond basic 
linear algebra, none of which are included below. 

Problems of the classical approach. As much as we all love the 
calculus, there have been limits to our applications in both pure and 
applied mathematics coming from the definitions which arose during 
the "rigorization" period of calculus. Leibniz had searched for an 'alge- 
bra of geometry', but this was not available until Grassmann realized 
the importance of /c-vectors in his seminal 1844 paper. But his ideas 
were largely ignored until Gibbs and Clifford began to appreciate them 
in 1877. Cartan used the Grassmann algebra to develop the exterior 
calculus. Meanwhile, Hamilton and others were debating the rigoriza- 
tion of multivariable calculus. It was not even clear what was meant 
by "space" in the late 19th century. The notions of div, grad and 
curl of Gibbs and Heaviside, were finally settled upon by the bulk of 
the community, and there have been two separate, and often compet- 
ing, approaches ever since. But the Cartan theory, rather than the 
coordinate theory we teach our freshmen, has given mathematicians 
and physicists a clearer vision to guide great leaps of thought, and 
this theory has led to much of the mathematics behind the prizewin- 
ning discoveries of mathematics and mathematical physics. The group 
of those who now understand the importance of the Cartan theory is 
growing, as evidenced by the many books and papers which now start 
with this theory as their basis. 

The reader might well ask, what are some of the problems of coordi- 
nate calculus? Does this Cartan theory have limitations? Why do we 
need something new? 

First of all, the coordinate theory requires multiple levels of limits 
upon limits to be able to understand interesting applications in man- 
ifolds or curved space that take years of training to understand. We 
become proud virtuosos of our coordinate techniques. But these meth- 
ods are a barrier to others who do not have the time or patience to 
learn them. And they also form a barrier to those of us who rely on 
them for they may cloud our vision with their complexity. After we 
become experts in Euclidean space with our div, grad and curl opera- 
tors, we then move into curved space. Everything works quite well in 
restricted settings of smooth Riemannian manifolds and submanifolds. 
We can manipulate the limits, know when to change their orders, and 
do term by term integration with Fourier series and wavelets, etc. But 
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when we try to reduce our assumptions and work with less smooth 
domains, permitting piecewise linear, corners, Lipschitz conditions, or 
introduce singularities, the coordinate theory becomes more and more 
intractable, and soon begins to break down. We are forced to consider 
too much information that is not necessary, and need more and more 
assumptions, and limits, to get anywhere. The classical calculus can- 
not treat everywhere nonsmooth domains as there are be no tangent 
spaces to work with. The Cartan theory cannot handle important bi- 
layer structures such as soap films because the boundary operator intro- 
duces extraneous terms along branched curves. The standard Cartan 
proof to Stokes' theorem relies on boundaries matching and cancelling 
with opposite orientation where they meet so local connectivity and 
locally finite mass are required. Discrete theories, so important today, 
also have problems. It can be surprisingly difficult to be certain when 
supposed approximations actually do converge to something meaning- 
ful. Bootstrapping techniques still prevail. Cochains in discrete or 
combinatorial theories usually fail to satisfy a basic property such as 
commutativity or associativity of wedge product, or existence of a sat- 
isfactory Hodge star operator. Simplicial complexes have appealing 
simplicity, but it is surprisingly difficult to model vector fields with 
them. There is a great deal of information in a simplicial complex 
that is not needed for calculus. There are corners, matching bound- 
aries of simplices, ratios of length to area, and so forth, that must be 
considered. 

As we try to apply calculus to important microscopic relations of 
physics or biology, the classical and discrete theories fail us. Both 
the Cartan model and the coordinate model assume a rigid, locally 
Euclidean continuum. This forces the calculus to completely break 
down in quantum mechanics. There is simple mathematics going on 
that cannot be seen with an assumption of a local Euclidean structure 
with it its "infinitesimal connectivity". 

Geometrization of Dirac delta functions. In this paper we present 
a new approach to calculus in which more efficient choices of limits are 
taken at key points of the development, greatly reducing the number 
of limits needed for the full theory. This work is motivated by Grass- 
mann's algebra of /c-vectors 

A = ©feA fc 

in a vector space V, but takes his ideas further. Instead of using k- 
vectors which have ambiguous geometric meaning, we do something 
slightly different to recover geometric meaning. We shrink a k-ce\\ to 
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the origin, say, renormalizing its mass at each stage so it remains con- 
stant. The question is "what do you get in the limit?" Intuitively, these 
geometrical "infinitesimals" are like "mass points", or geometric ver- 
sions of Dirac delta points. In this paper, we show that the limit exists 
in a normed space J\fg° and depends only on the fc-direction and mass 
of the k-ce\\. We call the limit a (0, k)-pole at the origin. The space of 
(0, fc)-poles at the origin is isomorphic to the space of exterior fc-vectors 
Afc. The boundary operator d : J\f£° — ► is continuous, and thus 

determines the boundary of a (0, fc)-pole as a sum of (1, k — l)-poles, a 
geometrization of Dirac dipoles of dimension k — 1. A geometric direc- 
tional derivative operator V u : J\f£° — > J\f%° is defined on (1, fc)-poles, 
leading to geometrizations of quadrupoles of abritrary order. We call 
V u a prederivative. Its dual is the directional derivative of forms. The 
formal treatment begins with the use of the Koszul complex X(V) in 
£0 Our geometrical view is supported and encapsulated by a matching 
algebraic construction of the Koszul complex which keeps track of all 
the implicit tensor algebra. For example, the Koszul complex shows 
that the boundary of a (0, fc)-pole is well defined, and leads to a nice 
proof that d 2 = 0. 

The real continuum beyond numbers. Arising from this theory 
are new models for the real continuum. Morris Hirsch wrote, ' on De- 
cember 13, 2003. 

A basic philosophical problem has been to make sense 
of "continuum", as in the space of real numbers, with- 
out introducing numbers. Weyl wrote, "The introduc- 
tion of numbers as coordinates ... is an act of violence". 
Poincare wrote about the "physical continuum" of our 
intuition, as opposed to the mathematical continuum. 
Whitehead (the philosopher) based our use of real num- 
bers on our intuition of time intervals and spatial regions. 
The Greeks tried, but didn't get very far in doing geom- 
etry without real numbers. But no one, least of all the 
Intuitionists, has come up with even a slightly satisfac- 
tory replacement for basing the continuum on the real 
number system, or basing the real numbers on Dedekind 
cuts, completion of the rationals, or some equivalent con- 
struction. 

We propose chainlets as more flexible models for the real continuum, 
since our "infinitesimal" (0, l)-poles can be limits of intervals, or equally 
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well intervals with countably many smaller intervals removed such as 
Cantor sets. The topology of limiting approximations of (0, l)-poles 
does not matter, contrasted with 1-dimensional tangent spaces which 
must have the structure of Euclidean space. Moreover, we add layers 
of possibilities to the continuum model, by allowing higher order poles 
at each point, a sum of a geometric monopole, dipole, quadrupole, etc., 
a kind of "jet" of geometry. Philosophers have said "there is only now" 
in human experience, but this continuum model gives a model of time 
which includes in each moment a jet of time, something we might sense 
as a Gestalt experience when time seems nearly to come to a halt, or 
the opposite when time seems to rush by. 

Our methods are different from those of nonstandard analysis which 
also defines "infinitesimals," but where there is great effort to mimic the 
structure of the reals as closely as possible. For the theory of hyperreals 
of nonstandard analysis to work, there is an ordering required, as well 
as inverses, transfer principles, etc. We only use local Euclidean space 
as an affine space to support our chainlets. We make use of the affine 
space structures of subtraction x — y and existence of norms \x — y\. 

Unification of viewpoints. Mathematicians take many viewpoints 
in the pure theory and its applications, be they smooth manifolds, 
Lipschitz structures, polyhedra, fractals, finite elements, soap films, 
measures, numerical methods, mathematical physics, etc. The choice 
sets the stage and determines our audience and our methods. All of 
these viewpoints can potentially be unified in chainlet theory. 

Three basic theorems for chainlet domains J lead to much of the 
classical theory, forming a "tripod" of calculus. 

(i) f 9J w = Jj dw (Stokes' theorem) 

(ii) JfjW = Jj f*w (Change of variables) 

(iii) f ±J w = Jj*w (Star theorem) 

Each of these is optimal, each has a one line proof, after the initial 
definitions and basic continuity results. (See §21 and |H4J.) Only the 
Star theorem requires a metric. With a metric, (i) and (ii) imply a 
general divergence theorem 




and a general curl theorem 
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In the standard approach, one proves the Fundamental Theorem of Cal- 
culus first, and then eventually builds up to a general Stokes' theorem 
for manifolds. We start with the three results above, as first results. 

In £1111 we provide a general method for testing whether a dual tri- 
angulation will converge to the Hodge star smooth continuum as the 
mesh size tends to zero. For example, this method, along with our 
change of variables result 15.91 and 15.61 can be used to show that the 
circumcentric dual introduced by by Marsden, et al, |D, H, L, M| does 
converge to the smooth continuum. (See fc llll for a general method.) 
The work in this paper, especially the continuity theorems for opera- 
tors, products and relations in §21 should fill in the missing ingredients 
needed to show which of the operators, products and relations in the 
above works converge to the smooth continuum. 

Besides the discrete theory, the author has developed two other ex- 
tensions of calculus: 

(i) Bilayer calculus with applications to the calculus of vari- 
ations including Plateau's problem (soap bubbles) (See |H8j . 

(ii) Calculus on nonsmooth domains (e.g., fractals) (See [HN1 , 

|HN2| . [TTTj . [H2], JH2, jHHO 

Differential complexes. Let Q r k denote the space of differential k- 
forms of class C r on a manifold M n . Since the exterior derivative d 
satisfies d 2 = it follows that 

ii fc+r "fc+i «fe "fc-l ii Q 

is a chain complex for each k + r = c > 1. Versions of this com- 
plex provide the basis to much of analysis and topology. Differential 
forms provide coordinate free integrands, leading to the Cartan exterior 
calculus De Rham cohomology theory is based on this complex. 
Classical operators such as Hodge star and pullback are closed in this 
complex, leading to broad applications. 

Mathematicians have sought a matching differential covariant com- 
plex for domains, with d playing the role of d and for which the complex 
is closed under smooth mappings. Poincare introduced the simplicial 
complex with its boundary operator mapping a simplicial complex of 
dimension k into one of dimension k — 1. Many papers continue to 
be written from the viewpoint of the simplicial complex, especially 
in numerical analysis. Although much progress has been made, (e.g. 
|D, H, L, M| ), the simplicial complex brings with it inherent problems. 
The smooth pushforward operator is not closed in the simplicial com- 
plex, making iterative methods problematic. Difficulties include the 
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commutative cochain problem, lack of associativity of cochains, lack of 
natural definitions of vector fields, a well-behaved discrete Hodge star 
operator, and convergence of operators and relations to the smooth 
continuum. We know of no other discrete theory, apart from that 
presented in this paper, which solves all of these problems simultane- 
ously. For example, wedge product is defined for our discrete cochains 
as restriction of the wedge product for forms at discrete points and is 
commutative, by definition. 

Beyond the needs of the discrete community, though, pure mathe- 
maticians have sought a geometrically based, covariant differential com- 
plex for which operators of analysis act continuously on the complex. 
Of especial interest are operators dual to operators on differential forms 
such as Hodge star *, and those which commute with the pushforward 
operator, for these lead to well defined integral relations in manifolds. 
An important goal has been to extend the Gauss divergence theorem 
to nonsmooth domains. (See Theorem 15.111 ) 

Whitney jW] introduced the vector space of polyhedral chains which 
does form a chain complex satisfying d 2 = 0. The Banach space of 
his sharp norm is not a chain complex because the boundary operator 
on polyhedral chains is not continuous. His flat norm does yield a 
chain complex with a well defined boundary operator. This has led to 
applications in geometric measure theory jFj. However, the Hodge star 
operator is not closed in the flat covariant complex. The flat norm has 
no divergence theorem, as can be seen by the example in Whitney of a 
flat form u in K 2 without flat components (jW], p. 270). 



The net divergence of u is zero, yet the net flux across the boundary of 
a square with diagonal x = —y is nonzero. The problem resides with 
the Hodge star operator which is not continuous in the flat norm*. 

Schwarz's distributions and de Rham's currents deRj solved the co- 
variant chain complex problem almost too perfectly. Currents T£ are 
defined to be continuous linear functionals on differential fc-forms of 
class C r with compact support. Bounded operators on currents are 
defined by duality with bounded operators on forms. For example, the 
boundary of a current T is defined by dT(u) = T(du). Stokes' theorem 
holds by definition. Since d 2 = 0, it follows that d 2 = 0, yielding the 



*More than one research group has not noticed this example, and tried to use 
the flat norm to develop calculus over fractal boundaries. 




dx + dy, x > —y 
0, x < —y 
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current complex 

T° A JL Tf-X & q-r 9 n-r+X 9 d rr-k+r 

2 k+r —>•••—> -l k+ i —> J- k ~~ * ^fc-l ~~ *■ ' ' ' ~^ 2 

Questions of regularity of currents can be a serious issue. For ex- 
ample, does the boundary of a current correspond to the boundary 
defined geometrically? Is a current solution smooth? What is its sup- 
port? These problems can be difficult to solve. In a way, spaces of 
currents are too big, and we look for proper subspaces that are sim- 
pler to work with. Examples include the normal and integral currents 
of Federer and Fleming [FFJ which have been important in geometric 
measure theory. 

We identify proper subspaces of currents J k <zT k r , called k- chainlets 
of degree r, that give us a natural covariant chainlet complex for do- 
mains that is algebraically closed under the basic operators of calculus. 

7 j> 9 rrr-X 9 rrr & rrr+X j> 9 ^fc+r 

While currents are dual to forms, chainlets are a predual to forms. 
More precisely, let B k denote the space of fc-forms, of class C r ~ 1+L%p , 
with a bound on each of the derivatives of order s, < s < r. (We 
do not require that the ambient space be compact.) Then the space 
of bounded linear functionals on J£ is precisely B k . Chainlets are not 
reflexive and thus J~ k C (J k )** = C T k r . An example of a current 

that is not a chainlet in V — R n is R n itself. However, even for compact 
ambient manifolds, chainlets form a proper subspace of currents since 
the space of chainlets is normed and the space of currents is not. 

Each differential fc-form u G B k is represented by a fc-chainlet in 
the sense that jj rj — f u A*n for all rj G B r k . (The LHS is the chainlet 
integral of £0 The RHS is the Riemann integral.) Therefore, B k is 
naturally immersed in J k . (See |H6j . |H7j .) 

All three of these spaces are Banach spaces for each < r < oo, and 
the inclusions are strict. Each is dense in the larger spaces. Define 
the direct limit = {J<7 k an d the inverse limit £>£° — [\B r k . The 
space £>jjf is a Frechet space, while J£° is a normed space. It is a direct 
limit of Banach spaces and contains a dense inner product space, while 
currents 7^°° are a topological vector space. 

The chainlet complex contains a dense subcomplex of discrete chains 
Vf.{V), called polypolar chains with each chain supported in finitely 
many points. 

-pO 9 d -pr-l 9 jj r d w+l 9 d -pr+1 

' k+r k+ i I k / k _i ■ ■ • r o 
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Wedge product and inner product are both defined on this discrete 
subcomplex. But there is no continuity of either of these products in the 
chainlet norms. The chainlet complex has neither a wedge product nor 
an inner product since it includes polyhedral chains and L 1 functions. 
(Polyhedral chains have no wedge product and L 1 functions have no 
inner product.) But the space of chainlets contains a dense, inner 
product space - the algebra of monopolar A;-chains. 
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1. The Koszul complex 

1.1. Exterior algebra. Let V be a vector space over a field F = Ror 
C. Let T(V) denote its tensor algebra. Let / be the two sided ideal of 
T(V), generated by all elements of the form v <g> v,v G V and define 
A(V) as the quotient 

A(V) = T(V)/I. 

Use the symbol A for multiplication, or exterior product in A(V). 
Since vAv = we deduce v Aw = —wAu and v±A- • ■ Av k = whenever 
W = (wi, . . . ,w k ) are linearly dependent in V. Elements of the form 
Vi A • • • Av k are called simple k-vectors. We sometimes denote these 
by 

a k (W) = tui A • • • A w k 
where W is the list (w±, . . . ,w k ). For example, a2(v,w) — v A w. If 
Wi = (twi.i, • • • , w lik ) and W 2 = (w 2 ,i, w 2 ,i), let 

(Wi, W 2 ) = (w hl , . . . , w ltk , w 2 ,i, • • • , w 2 ,i). 

Thus akiWx) A ai (W 2 ) = a k+l {W u W 2 ). 

The subspace of A(V) generated by all simple k- vectors is known as 
the k-th exterior power of V and is denoted by A k (V). The exterior 
algebra can be written as the direct sum of each of the k-th powers: 

oo 

A(y) = 0A fe y 

fc=0 

noting that Aq(V) = F and Ai(V) = V. The exterior product of a j- 
vector and an A;-vector is a (j + /c)-vector. Thus, the exterior algebra 
forms a graded algebra where the grade is given by k. 

We may replace v 2 = v <g>v with any quadratic form Q(v) and obtain 
geometric differentials of the Clifford algebra. This extension will be 
treated treated in a sequel. 

1.2. Symmetric algebra. The symmetric algebra of V is defined 
as follows. Let J be the two-sided ideal of T(V) generated by all 
elements of the form v <g) u — u <g) v. Define S(V) = T(V) / J. The 
symmetric product of u, v G V is denoted uv = vu. For j > 1, 
denote the j-th symmetric product by 

V 3 V = u\u 2 ■ ■ ■ uj 

where U is the list (ui, . . . , Uj). The subspace of S(V) generated by 
the j-th symmetric products is denoted S^(V). There is a direct sum 
decomposition of S(V) as a graded algebra into summands S(V) = 
®S j (V) where S°(V) = F and S\V) = V. The space S j (V) is the j-th 
symmetric power oiV . Elements of S j (V) of the form U\U 2 ■•• Uj,j > 
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1, are denoted by V{j where U = (u±, 112, ■ ■ ■ , Uj). We observe that V{j is 
independent of the order of the vectors of U and permits duplications, 
whereas a k (W) depends on the order of W . 
Define the associative, unital algebra 

X(V) = ® J , k S>(V)®A k (V). 

Denote V 3 v a k = V 3 V <g> a k , where G S j (V) and a k G A k (V). The 
integer k is the geometric dimension of V 3 u a(V k ) , and j is its order. 
The product ■ in X(V) is defined for simple vectors a(Wi),a(Wi) by 

<a(Wi) • Vg,a(W 2 ) = V^ a) o(Wi, ^ 2 )/2^ 2 , 

f° r ji) J2 > 1) and extended by linearity. For j\ — or j 2 = 0, the 
product is merely scalar multiplication. The operator V\j : X(V) — > 
X(V) is defined by 

V^V&a) = Vf+^a. 

1.3. Mass and direction of simple A>vectors. We choose a pre- 
ferred basis (ei, ... , e n ) of V and use it to define an inner product on 
V: < e iy ej >= 5^. Define mass of a simple /c-vector by 

M(a(w±, . . . , tUfe)) = det(< Wi, Wj >). 

If a different basis is chosen, the two masses will be proportional. The 
resulting normed spaces will be the same, and the theory will be iden- 
tical. The orientation of the subspace W of V spanning the list 
(wi, . . . ,w k ) is the orientation of the list (wi, . . . ,w k ). (Two linearly 
independent lists have the same orientation if and only if the determi- 
nant of the change of basis matrix is positive.) The k-direction of a 
simple /c-vector a(wi, . . . , w k ) is the oriented /c-dimensional subspace 
spanning the list (w±, . . . ,w k ). The fc-direction of a simple k- vector is 
independent of the preferred basis. 

Lemma 1.1. Two simple k -vectors a and (5 are equal if and only if 
their masses and directions are the same. 

The proof follows directly from the multilinearity properties of tensor 
product. 

1.4. The Banach algebra X(V). Define 

|| V^-aHj = |«i| • • • \uj\M(a) 
where U = (ui, . . . ,Uj). A basis of V generates a basis of ® A*,. Let 
A 7 = J2i a i^u a ii wr hten in terms of this basis, and define 

ii^l- = Ewil v Wli- 



12 



J. HARRISON DEPARTMENT OF MATHEMATICS U.C. BERKELEY 



Finally, for A = £\ A? G S(V), define 

\\A\\=J2W AJ h- 

j 

Proposition 1.2. || • || is a norm on the algebraX(V) satisfying 

\\A-B\\<\\A\\\\B\\. 

Proof. Supppose A = J2j=i ^ an d ll^ll > 0- We may assume that 
the terms of A^ = Yli a i^Ui a i are written in terms of a basis. By the 
triangle inequality, ||^4 J ||j > for some j. Hence || V^.a^Hj > for some 
i. By definition of the norm, and since mass is a norm, it follows that 
V^.tti 7^ 0. Since A* is written in terms of a basis, we know A 3 ^ 0, 
and thus A ^ since we are using a direct sum. The other properties 
of a norm follow easily from the definitions. 

The inequality is a consequence of the triangle inequality and 

HVj^i • VfaiWj = \\D$:$ t) ai A ael 
= \Ui\\U e \M(a t Aa e ) 
< \U t \\U t \M{ai)M{a,) 

= \\VuM\ii\\ V uM\ir 

□ 

1.5. The covariant complex. 

Prederivatives V„. Let u e V and V^a G -X"(V). Define the pred- 
erivative ( in the direction u) by 

V u : S j ®A k ^ S^ 1 <g) A fc 

by 

V U (V^) = Vg^a. 

In the next section we give a geometric interpretation of this operator. 
In a sequel, we show it is a derivation. 
We next define the boundary operator 

d : S j <g> A fe -> S^ 1 <g> A fc _i, fc > 1. 

Define 

^:=V,(1). 

There is a unique extension of d to -X"(V) making the boundary operator 
into a derivation. In particular, 

d(u Av) = du ■ v — u ■ dv. 
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In general, 

8 : S j <g) A k -> <g> A fc _! 

is defined recursively: 

Assume the boundary of a simple /c-vector a lias been defined, and 
v E V. Define 

d(a A v) = da ■ v + {-l) k a ■ dv e S 1 <g> A k . 
For j > 0, define 

0(v£a) = Vi(Sa). 
We obtain a linear mapping 

d : s-J'-^y) ® A fc (y) - s*(v) ® A fc -i(^) 

which is a derivation. It follows that 8 o 8 = 0. 

We conclude that is a unital, associative, bigraded differential 

Banach algebra. The product of X(V) is associative, bilinear, graded 
commutative, and the boundary operator is a derivation: If A, B e 
X(V), then 

8(A ■ B) = (OA) ■ B + (-l) dim W ■ (8B). 

For each < k + j = c, < k < n and j > 0, we have the bigraded, 
differential, covariant complex 

s°(v) ® A j+fe (y) ^ s 1 ^) ® A j+fc -i(t0 ^ • • • ^ ® a (v). 

We next see how the boundary operator relates to the prederivative 
operator in the case of a simple fc-vector. 

Lemma 1.3. If a is a simple k-vector, then da is the sum of k pred- 
erivatives of simple (k — l)-vectors. 

Proof. The proof proceeds by induction on k. The result holds by 
definition for k — 1. Assume it holds for simple (k — I)— vectors. If a 
is a simple k- vector then a — (3 A v where (3 is a simple (k — l)-vector, 
then da = 8(3 -v + (— l) k ~ l (3- dv . By induction, 8(3 is the sum of (k— I) 
prederivatives of simple (k — l)-vectors. But (3 ■ dv = V v (3 is also the 
prederivative of a simple (k — l)-vector. □ 

1.6. The contravariant Koszul complex. Define the exterior de- 
rivative operator 

d : S*(V) ® A k (V*) - &-\V) ® A k+1 (V*) 

as follows: d(i> ® 1) = 1 ® u. Much as in the preceding section, d has a 
unique extension so it is a graded derivation on the product. 
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For each < k + j = c, < k < n and j > 0, we have the contravari- 
ant Koszul complex 

s°{v) ® A j+k (v*) s j+k -\v) ® Ai(y*) 4- ^ +fc (vo ® a (v*) 

Theorem 1.4. ® A fc (V*) c (S^V) <g> A fe (y))*. 

Even though A k (V*) (A fe (V))* and ^ the tensor 

products of dual spaces is not isomorphic to the dual of tensor products. 
The exterior derivative that we have defined as dual to the boundary 
operator is an operator on the space (S^V) ® A k (V))*, and thus is 
also an operator on ^'(V) <8> A k (V*). On the other hand, the exterior 
derivative defined in the contravariant Koszul complex 5 J '(V) ® A k (V*) 
by "moving covectors from the symmetric side to the antisymmetric 
side" is not naturally extendable to (S j (V) (g> A k (V))*. 

The algebra of polynomials is isomorphic to the symmetric algebra 
and can be useful as coefficients of fc-covectors in S^(V) <S> Ak(V*). 
Again, this does not have natural extension into (S^iV) <8) A k (V))*. 
For (S^(V) ® A k (V))* it is more natural to use differential forms with 
coefficient functions with bounded derivatives, such as trigonometric 
functions or smooth functions with compact support. This touches on 
a philosophical difference between Taylor's theorem approximations by 
polynomials, vs approximations with Fourier series or wavelets with 
their superior convergence rates. 

In summary, the covariant complex X(V) leads to a more general 
theory than the contravariant complex X(V*). It is simpler to work 
with conceptually as it is based upon physical concepts. ^ The basic 
operators of pushforward, boundary, and perp are more natural than 
the duals of pullback, exterior derivative, and Hodge star. 

2. The polypolar chain complex 

2.1. The algebra of polypolar chains. Consider the product space 
V x X(V). A j-polar k-vector in V with support at p is a pair 
(p; V^a) G V x A 3 k (V). We form the vector space V J k (V) of j-polar 
k-chains of formal sums A = Yli=i(Pi'> ^u a i) subject to the relation 
ip'i ^u a ) = (Pi — ^u@) if a an d P have the same mass and fc-direction, 
but with opposite orientation. Furthermore, (p; V^a) + (p; V^,o/) = 
(p; W j v a + V{j,a'), and £(p; V^a) = (p; tW j v a),t e F. Define V k = ®jV j k 

^There is a foundational difference between the approach of |T] and our own. By 
building upon the prederivative operator V u of the Koszul chain complex, rather 
than the exterior derivative operator of the Koszul cochain complex, we generate a 
discrete theory with a geometrical basis rather than a basis of differential forms. 
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and V = ®kPk- For j = 0,1,2, we call j-polar /c-chains monopolar, 
dipolar and quadrupolar A;-chains, respectively. We are especially in- 
terested monopolar fc-chains as they are dense in spaces of chainlets. 
For simplicity, consider monopolar 1-vectors. The space Vi(V) differs 
from the affine space of V in the following important way. The space 
of monopolar 1-chains is a vector space, equipped with the translation 
operator T v {p; V u a) = (p + v; V u a). The affine space associated to a 
vector space V is not a vector space, as there is no addition, only sub- 
traction p — q. A monopolar k- vector (p;a). separates the two basic 
roles of a vector into its two parts. In the first slot of (p; a) , we think 
of the vector p e V as a point. It is merely the support, or location, of 
the fc-vector a in the second slot. The pair becomes a "smart point" 
which contains much more information than just p alone. We make use 
of two projections: 

supp : V J k (V) — > V, 
mapping P = J2(pi] V^.ctj) to the point set supp(P) = Upi and 

Vec{ : Vi(V) Ai(V) 

mapping P = X)(p»I V^.aO to V^.a,. 
Wedge product on V . Define 



(PI a) A (q; (3) 



0, p^q 
(p;aA/3), p = q. 



For A = J2( Pi ; AA and B = ^(qf Bj) G V x X(\/), define 

where = qj. This reduces to wedge product of A;- vectors if A and B 
have order zero. For A = ^(puAA, define the norm 

H| = £11411- 

Then 

||AA5||<|LA||||5||, 

making the space of monopolar chains V into an associative, unital, 
Banach algebra. 

Define the 1-difference monopolar k -vector by 

A l(p; «) = T u(p; a) - (p;a) 

and the j -difference monopolar k-vector recursively by 

Ai,(p;a)=Ai(Aj7 1 (p;a)). 
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where U = (u,U'). Define a;(A*(p; a)) = uj{T u {p;a)) — uj(p;a) and 
extend recursively to uj(AIj(p; a)) and by linearity to uj(P), P E V{. 

We next see how to develop calculus at a point by way of polyvectors. 
Define 

||A°(p;«)||o = M(a), 
and, for j > 1, define ||A^(p; a)\\j = \U\M(a) where \U\ = \u±\ ■ • • \uj\, U = 

(«!,••• ,Uj). 

Define T> 3 k (V) to be the subspace of V®(V) generated by j-difference 
/c-chains 



2.2. Differential forms on monopolar chains. A differential k- 
form uj is defined to be a linear functional on the space of monopolar 
/c-chains Vk(V). The support of a /c-form is a closed set K defined as 
follows: A point x is in the complement of K if there exists a neigh- 
borhood U of x and missing K, such that u>(x; a) = for all pairs 
(x;a),x eU,a e A k (V). 
Denote the operator norm 

\u;\h = \\u;\\ = sup {\lu(x; a)\ : M(a) = 1} . 

We say that uj is bounded measurable if \ujy° < oo. The space of all 
bounded measurable forms is denoted Vk{V)* 
Let uj e V k (Vy and j > 0. Define 



= sup{\uj(Af(x;a))\ : || Af(x; a) \\ 3 = 1} . 
Define the norm 

{uj^ = max{\\uj\\ , . . . , \\oj\\ r } 

Proposition 2.1. \uj\^ r is a norm on the subspace of bounded Lipschitz 
forms with \uj\^ r < oo. 

The completion of the space of bounded measurable /c-forms with 
Icj]^ < oo is denoted B 7 k (V)} Examples include sin(2x)dx, f(x)dx 
where / has compact support. Nonexamples include xdx since ||xdx||o = 
oo. 



*In coordinates, this space is equivalent to the space of fc-forms with coefficients 
of class C r ~ 1+Llp , i.e., the Lipschitz constant of each derivative of order (j — 1) is 
uniformly bounded in V. 
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2.3. Wedge product of forms. Given lo G B r k and rj G B s m , define 



LuAr)(p; a(v 1: . . .,v k+rn )) 

= 5^ u ^ • • • ' v <y{k))r]{a{p] v a{k+1) , . . . , v a{k+m) ). 

<T&S k 

Proposition 2.2. Suppose u e B r k ,r/ e B s m . Then 

\uj A r]\^ +s < \oj\^ r \r)\^ s . 

Proof. We prove that 

||W A ^llt+ti < IMItlMU 

for all < t < r, < w < s. The LHS is the supremum of terms of the 
form ( ^?li Ay °"° +m \ We may assume wlog that afc +m = a(e±, . . . , e k+m ) 

W^ZJ a k+m,\\t + u 

where (ei, . . . , e k+m ) is orthonormal. Then M(a fc+m ) = M(a fc )M(a; m ) 
for any permutation of the vectors of a k+m . Split the list JJ t+u into two 
lists U t+U = (Ul,U%). 

The result now follows since wedge product is averaged over all such 
permutations. □ 

2.4. Norms on P%(V). Denote the dual norm on V k (V) by 

|P|^ = su P ^l, P eV° k (V). 

By duality, \P\^ r is a norm on the space V k (V). 
For D j G W k (V), define 

II ni II U ^ Dj) 
\\D J \\j = sup — — - — . 

\M\j 

The next lemma follows directly from the definitions. 

Lemma 2.3. |£tf|k < \\D^\\j. 

Theorem 2.4. IfPeV% and r>0, then 

\P\^ = inf^jDi, : P = YT j=Q Di,Di G V{}. 
Proof. Suppose P = . By the triangle inequality 

\pf r < ^ \ Dj \^ r < w^Wi- 
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On the other hand, let e > 0. There exists P = ^2 such that 
RHS > £ \\D3\lj -e. Then 

\Py r = SUp - — = SUp 



^ I ih , .. M M || 

ogBI \wr we i%m&x{\\u\\o,...,\\u\\ r } 

f io(D°) uj(D t 
< sup <^ + ■■■ + 



ojgbi L IMIo 1Mb 

r 

< IPii < RHS + e. 

j=0 

Since this holds for all £ > 0, the result follows. □ 
2.5. Prederivative operator V w . 

Lemma 2.5. If lu e Bl, p <E V, and a e A fc (V), t/ien 



D u uj(p; g) = lim u; I p 



exists anc? is unique. 

The proof to this appears in [H6] . If a; is only Lipschitz, Radamacher's 
theorem shows this limit exists a.e. If uj G i3 1+£ , the limit exists and is 
continuous. 

Define the directional derivative 

Afa(p; a) 

and 

w(p; V u a) = D u u(p; a). 
Lemma 2.6. D u is bilinear in u and a and satisfies 

\D u u\^- x < \u\\lo\^. 

Proof. It suffices to show 

||£) u a;||j_i < \u\ 

But 

n „ (P^XA^fog)) ^(A^faaQK 

MAl ,_i = sup —-, — = sup lim t — < \u\ \\uj L- 

' '|A?r x (p;«)lli-i ^° i||A| Ui[/) (p;«)||- 

□ 

Lemma 2.7. |V U P|^ < lullPl^- 1 . 
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Proof. By Lemma 12.61 

V K P tr = sup ," ; = sup " , ' < sup 1 , ,,' 1 < Lu LP 

\u\* T \uir r \uir r 

□ 

2.6. Exterior derivative d. By Lemma fl. 31 the boundary of a simple 
k- vector a is the sum of k prederivatives of simple (k — l)-vectors 
Pi. That is, da = Y^^uiPi- It is often convenient to assume that 
the directions of translation are orthogonal to the This follows if 
a = a(vi, . . . ,Vk) where the list (vi, . . . , Vj.) is orthogonal. But a always 
has such a representative by Lemma fl.ll 
Define the exterior derivative on forms by 

du(p; a) = uj(p; da). 

Lemma 2.8. Let u G B\ and <j) e B\. Then 

(i) d(f)(p; u) = V u 0(p) 

(ii) du(p;a) = EV»,w(j);A). 

Proof. Part (i) follows since 

d(f>(p; u) = (f>(p; du) = <f>{p; V u (l)) = D u (p(p). 
For Part (ii), we have 

du(p; a) = u>(p; da) = u(p; £V Ui $) = D Ui u(p; Pi). 

□ 

Proposition 2.9. (i) d(u + r]) = dw + drj 

(ii) d{oj At]) = did A r? + (-l) k LU A dr/ 

(iii) d{du) = 

(iv) d(p = D e .(j)dei if (ei, . . . , e n ) is orthonormal. 

Proof. Part (i) follows from linearity of boundary operator. Parts (ii) 
is a consequence of Lemma 12.81 and the definition of wedge product. 
(See (FT], for example.) (iii) follows since d o d = 0. 
Proof of (iv): 

□ 

Proposition 2.10. If u e B r k , then 

\dw\^-i < {k + l)\uj\^. 
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Proof. Let a be a simple (k + l)-vector. Suppose da = J2 V Ui A where 
Ui is orthogonal to the /c-direction of Therefore, 

\\du\\j-i = SU P I, A3 --i n — = SU P NAi -i n — 

ll A l/ ll A [/ a \\j-l 

^ E^(AJ Mi[/) A) 

< sup 



< (* + l)IM 

= (* + i)IM 



«i||C/|M(A 



|t/|M(AA^ 



□ 



We may now extend the definition of forms recursively so they are 
defined on each V r k . Assume that forms of class B k +1 have been defined 
on V J k ,j < r. If uo E B' k +1 , then D u uo G B\. Define 

^ r u,u a ) = AMp; V^ -1 ")- 

We may now define 

du(p; V^a) = u(p; V^da). 

The boundary operator extends by linearity to monopolar chains P. 
Stokes' theorem for monopolar chains follows 

doo(P) = oo(dP). 

Proposition 2.11. If P e V r k , then \dP\^+ l < k\P\^. 

Proof. By Proposition 12.101 

nnl . , u(dP) du{P) , \dcuMp\br 

\q p Wi = sup \ > = sup , , ; < fcsup 1 , ' 1 = klPr. 

□ 

By linearity dod = 0, giving us a chain complex for each k + r = c : 

-pO d 9 -pr-l 9 -nr J* -pr+1 9 9 -pfc+r 

' k+r k+ i I k / k _i rg 



chainlets 
3. Smooth Pushforward 
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We first give a coordinate free definition of the r-norm on smooth 
mappings, and then relate this to a more familiar norm relying on 
coordinates. 

Definition 3.0.1. The r-norm of a mapping f : V — > M n is defined by 

\f*Aij(x;a)^ 



[j] = SU P 



and 



\f\[r] = max{||/|| [0 ],...,||/|| [r ]} 
We say f e B r if\f\ [r] < oo. 

We remark that ||/||[i] < C is equivalent to assuming that / satisfies 
a Lipschitz condition. 

\f(x + u)-f(x) 



< C 



\u\ 



Let Qk(0) denote the n—cube centered at the origin with edge length 
IK. 

Lemma 3.1. Let f = . . . , f n ) : V — > Qk(0) in standard coordinates 
ofR n and r > 0. Then for each 1 < % < n, 

\fit<K\f\ r . 
Proof. It suffices to show each 

UWi < ll/llbi 

Let 



Ui(x; 1) 



\xA < K 



K, \ Xi \ > K. 



Then 
But 



l^illo < K, \\cu\\ s < 1 



\wi\v < K. 



l/i(Ai(x;l))| 



fi 



A'(s;l) 



Therefore, 



f»A 3 u (x;l) 

<\u i f>\f m A>,(x;l)f 
<K\f*Aij(x;l)^ 

ll/i|li<^ll/llw =► l/il 11, < ^l/lr- 



□ 
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Lemma 3.2. Let f e B 2 ,u &V,x<eV. Then each directional derivative 
D u fi(x; 1) exists and is a 0-form of class B 1 . 



Proof. 



DJi(x; 1) = fi(x; V u (l)) = Jim fi(A tu (x; l/t)) 



by the main lemma. □ 
For u = ej, then D u fi(x; 1) = 

Definition 3.2.1. Let f = (fx, . . . , f n ), x E V, and u £ V . Define the 
pushforward f*(x;u) by 

f*(x; u) = (f(x); D u fi(x), ■■■ , D u f n (x)) 

One way to interpret pushforwards is as follows. Let / : R m — > R n be 
a smooth mapping. Then we can think of f* x : R m — > K m as the best 
linear approximation of / at x. Suppose u e K m . Let <7 U = {x + t-u : 
< t < 1} be the 1-cell of (x;u). Then f*a u is a singular 1-cell with 
endpoint f(x). The tangent vector to f*cr u at /(x) is f*(x;u). Now 
/* x is linear and in coordinates takes the form of the classical Jacobian 
matrix 

••• §t^\ 

In order to extend to monopolar /c-chains, we define 
f*(x; tii A • • • A Wfc) = (/(x); A • • • A 

Define 

/.(a;; V u a) = (/(x); V/.„/,,a) 

From here we provide the following useful lemma that relates push- 
forward to all the other operators we have studied: 

Lemma 3.3. 

(i) /*V U = V/* u /* 

(ii) (fog)* = 9*° f* 

(iii) /*d = df* 

Note however that the pushforward and perp operators do not com- 
mute. A counterexample is the shear operator. 

Another property that is standard is the chain rule. This follows 
since it holds for linear mappings. 
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Theorem 3.4 (Chain Rule). Suppose f : R m — > M n and g : R n — > R p 
are smooth mappings. Then 

(f°g)* = f*g*- 

Proof. 

f*g # (x;v) = f*(g(x); g* x v) = (f(g(x)); f* g ( x )g* x v) 

= (/ ° 0(3); (/ ° = (/ ° g)*fa v ) 

□ 

Definition 3.4.1. Let f : R n R p be a smooth mapping and uj a 
k-form on R p . Define f*u on M n by 

f*u(x; a) = uf m (x; a) = uj(f(x); 
The change of variables theorem is immediate for monopolar chains. 
THEOREM 3.5. Change of variables: 



f*u= u 
J }*p 

The pullback behaves similarly to the pushforward. 

Lemma 3.6. 

(i) f*{uAr])=f*u J Af*r] 

(ii) {fogf = g*of* 
(hi) f*d = df* 

Proof. 1. follows since it holds for covectors. For 2. we apply the chain 
rule. 3. follows from Lemma f3. 31 □ 

Now we establish the continuity of the pushforward operator. 

Theorem 3.7. \f*w\* r < \f\[ r ]\^ r 

Proof. 



\ A u(x;®)\\j ' \\ A u(x;ai) 



< wf j sup - — y 

II A-? 



< \f\[r]\^ 

The result follows. □ 
Now we can establish the continuity of the pushforward. 
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COROLLARY 3.8. If P is a monopolar chain and f is a smooth mapping 
of class B r , then \f m Pfr < \f\[ r ]\P\* r 



Proof 



. KAP)j |/*g;(P)| 

\f*Pr = SUp —r- = SUp [—7 < SUp 

\uir r \U) r r 



\f*co 




\P\ 


tfr- 




Jj\ 


3r 



<sup "'" 1 1 =\f\ M \P\*- 

\Ur r 

□ 

4. The operator _L 

Let a be a simple A;- vector. Define 1 a to be the simple (n — k)- 
vector with the same mass as a and with (n — fc)-direction orthogonal 
to the fc-direction of a. Define 

JL (V tt a) = V„(± a). 

Extend by linearity to define 

J_. T>J -p-?' 
- 1 - ' fe ' n-k- 

(Sec [H5J for details of the perp operator in R n .) 
Define the unit volume n-form dV by 

dV(a) =J_ a 

for simple n- vectors a. By linearity f p dV =J_ P for P G V®. 

Define the unit volume n- vector by vol = a(ei, . . . , e n ) where (ei, . . . , £ 
is an orthonormal basis of V. 

Lemma 4.1. (i) _U_ a = (-l) fc ("- fc )a 

(ii) aAla = M(a) 2 vol 

(iii) J_=J_ 0. 

Proposition 4.2 (Continuity of the perp operator). If P e V r k , then 

Define d* = =_L <9 J_ and the geometric Laplace operator by 
□ = d o + o 9. Define the Hodge star operator by 

*uj{p; a) = u(p; J_ a). 

Lemma 4.3. 

I *U)\^ r = \u\^. 
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Proof. 

^(A^faa)) u(Ajj(p; JL a)) 
*w ,■ = sup ^ = sup % 

U ^\\j\\ A u(P^ a )\\j- 



< 



|AMp;a)||, 



□ 



Theorem 4.4 (Star theorem). IfueB r k and P G V r n _ k , then 



~kUJ — UJ. 
P J±P 



The Laplace operator on forms is given by A = d5+5d where 5 = *d*. 
Proposition 4.5. Au(P) = d(aP). 
Theorem 4.6 (Star theorem). 



-kUO = I LU. 
P J±P 

5. The chainlet complex 

Denote the completion of the space of the space V k (V) under the 
norm \P\^ r by J\f[{V) and call its elements r-natural k-chainlets in 
V. 

Compare this to the definition of currents T£ with norm 

\Tr r = sup , . 

Currents are dual to differential forms. That is, (BQ* = T k . Chainlets 
of Ml are a predual to differential forms in that (Ml)* = B k . The spaces 
are not reflexive since M k is separable. It contains a countable dense 
subspace whereas, B r k is not separable. (See Whitney [W].) Therefore, 

mi c wr = n- 

Let J G M k . Then there exists a sequence of monopolar chians 
Pi — > J in the r-norm. Since the norms are decreasing, it follows that 
Pi converges to a unique chainlet f rs (J) G Ml +S , s > 0. 

Lemma 5.1. TTie mapping 

f„ : Afc - A/^ 

is a homomorphism. Furthermore, f rr {J) = J for all J G Ml and 
fn = fst ° frs, for allr <s < t. 
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Define the direct limit M^{V) = U7V fc r (V). 
Theorem 5.2. A/^° is a vector space with norm 

\P\*°° = lim|P|^. 

Proof. Let P = Yli=o(P^ a i) ^ e a monopolar /c-chain that is not zero. 
Then (p ; ao) ^ 0, say. We may assume without loss of generality that 
the fc-direction L of ao is the subspace spanned by (ex,... ,ek). Let 
rj = dxx- ■■dx k . Now r)(P) = YlUoViPii «*) = Y.Uo'tiPu A) where A 
is the projection of a» onto L. Since $ is top dimensional, there exist 
scalars f(pi) such that $ = f{pi)e\/\- ■ -Aet- There exists a polynomial 
Q defined on L such that J2Qi.Pi)fi.Pi) > 0- Define = Qr/. Then 

c(p) = gr / (p) = ^gfe)/fe)>o. 

Now each derivative of u; is bounded above by the derivatives of Q. It 
follows that |P|^ > 0. 

The other properties of the norm are straightforward to verify. □ 

Every monopolar fc-cell (p; a) is naturally included in Af[ for each 
r. The pair is called a (0, k)-pole. This corresponds naturally to an 
electric monopole for k = 3. 

The inverse limit £>jf(V) = nBl(V) of spaces of differential forms is 
the space of smooth A;-cochainlets. These forms are infinitely smooth 
with a bound on each derivative of order r. 

Lemma 5.3. B™(V) is a Frechet space. 

Theorem 5.4. Chainlets form a proper subspace of currents. 

Proof. This follows since chainlets form a normed space and currents 
are just a topological vector space. □ 

Chainlets are a rich and interesting space. They are more than just 
a normed space. They are the direct limit of Banach spaces, and have 
an orthonormal basis. So they are "almost a Banach space". They 
contain a separable, dense inner product space, and are thus "almost 
an inner product space" . 

5.1. Isomorphisms of forms and cochains. By definition, B r k = 
(A/^)*. Therefore forms and cochains are canonically isomorphic in 
this category. The isomorphism respects the pullback operator since 
df* = f*d. (See [H2J for a much longer proof of the same result.) 

THEOREM 5.5. The spaces B r k and {Ml)* are naturally isomorphic. 

The de Rham isomorphism of cohomology will later be found as a 
corollary since singular chains are dense in chainlets in a manifold. 
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5.2. The part of a chainlet in an open set. Let U C R n be open. 
For a monopolar chain P = ^(puca), define P|_t/= J2(.Pij' a ij) where 
Pij are points in the support of P D U. (See [H6J for more details of this 
section, including a discussion of nonexceptional open sets.) Suppose 
J has finite mass. If U is nonexceptional w.r.t J then if Pj — > J with 
M(Pj) — > M(J) then Pj|_£^ is a Cauchy sequence. We denote the limit 
by J [17. 

5.3. Improper chainlets. A sequence of monopolar chains Pj is said 
to converge to an improper chainlet J if Pi\u converges to a 
chainlet, denoted J\_u, for each open set U in V. An improper chainlet 
if is the boundary of an improper chainlet J, if d( J[y) — J[du= K\_u for 
all nonexceptional U. We write dJ = K. An example of an improper 
chainlet is R n . Consider a binary mesh of E n with vertices Let Q 
denote the unit n-cube and define Pj = ^2 _rifc (j>j; Q). It is left to the 
reader to verify that Pj converges to the improper chainlet R n . 

5.4. Operators on chainlets. Operators on monopolar chains are 
continuous in the chainlet norms, and therefore extend to operators on 
chainlets. 

5.4.1. Prederivative operator V„. It J e Ml is a chainlet, choose Pj — > 
J in Af[. According to Theorem 12.71 we know V u Pj is a Cauchy se- 
quence. Define V U J = lim V u Pj. Therefore, V u : J\f£~ — > A/£ satisfies 

|V U J|^ < lullJl"- 1 . 

The other operators and relations described in Section 3 for monopo- 
lar chains extend in a similar fashion, including boundary, pushforward, 
prederivative and perp. 

5.4.2. Boundary dJ. Define dJ = lim<9Pj. Since the boundary oper- 
ator is continuous, it follows that d o d = on chainlets, giving us a 
chainlet complex 

K,(Y) 4 • . . 4 M£Xy) *> mxy) 4 M£1(V) 4 ■ ■ ■ 4 AT k+r (v). 

5.4.3. Integration over chainlet domains. Define the integral over chain- 
lets J G A/2 - (V) and for forms u is of class B r k as follows: 

I to = lim u(Pi) 
J J 

where Pj — > J in the chainlet norm. The limit exists for J e Af[ and 

ijj G Af[, < r < oo, since 

\u(P)\ < \u\^\P\^. 
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The support of a chainlet is defined by its complement, using the 
chainlet integral. A point p is in the complement of supp(J) if there 
exists a neighborhood U missing K such that JjU = for all forms 
supported in U . 

By taking limits, we deduce the following results. 

Theorem 5.6. If J £ Af[, < r < oo, and uj e B r k , then 

I u < \J\ kr \uj\^. 

The next result is a primitive form of Stokes' theorem as it relates 
prederivatives of chainlets to directional derivatives of forms. 

Theorem 5.7. For u e B r k +1 and J e A^ r , < r < oo, 

D u uj = uj. 

J JVuJ 

The following general form of Stokes' theorem extends the classical 
version and applies to all chainlets J, including soap films, polypolar 
chains, and fractals. 

Theorem 5.8 (Chainlet Stokes' theorem). For uj e B r k +1 and J E 



oo, 



/ uj — du. 

JdJ J J 



Proof. 

du = lim / duo = lim / uo = I 

J J J Pi JdP t Jd 



UJ. 

i, i 

□ 



Theorem 5.9 (Chainlet change of variables). For uj e B r k , J e j\f k , 
and f G B r , < r < oo, 



Proof. 



f*U = / UJ. 

J J ft J 



f*uo = lim / f*uj = lim / uj = uj. 

J JP t Jf*Pi Jf*J 



□ 



THEOREM 5.10 (Chainlet star theorem). 



f*u= f 
J j J i 



UJ. 

±J 
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Corollary 5.11 (Chainlet divergence theorem). 

d-k LU — 10. 
'J J±dJ 

Corollary 5.12 (Chainlet curl theorem). 



•kduO — I LU. 
I J Jd±J 

5.5. Four dense subsets of chainlets. 

5.5.1. Monopolar chains. Monopolar chains are dense in chainlets, of 
course. Monopolar chains V are not only a Banach algebra but they 
are also an inner product space, both structures not present in chain- 
lets. But the inner product gives us certain numerical advantages for 
chainlets, since an orthonormal basis for V is also an orthonormal basis 
of each J\f£, < r < oo. 

Inner product on monopolar k-chains. Define 

<p,g>= / dV. 

JPA*Q 



The norm \P\ 2 = V< P,P> = a/S M(di) 2 is an "L 2 norm" for 
monopolar chains. For an 11 LP norm", p > 1, define 

\p\ P = C£ M ^) p ) 1,p - 

One may extend the geometric product of Hestenes simply by taking 
the sum PQ =< P,Q > +P AQ, P,Q 6 V%. Then PQ is a well-defined 
element of the algebra of monopolar chains V. It is an open question 
whether this will lead to anything more than chainlets already present. 
A forthcoming extension of the theory to Clifford algebras, however, is 
clearly deeply important. 

5.5.2. Polyhedral chains. Polyhedral chains were invented by Whitney 
and were the basis to his Geometric Integration Theory jWj. We show 
that polyhedral chains are naturally included in the space of chainlets 
as a dense subset. 

Proposition 5.13. A k-cell a naturally corresponds to a chainlet. 

Proof. First consider a /c-cube Q. Divide it into subcubes Qi. Let Pi 
be the midpoint of Qi. Let (pf, ai) be the unique (0, fc)-pole supported 
in pi with mass the same as the mass of Qi and the same fc-direction 
as that of Q. The monopolar chain XXp«> a i) ls Cauchy in the natural 
norm. (See |H6] .) Call the limit Jq. Finally, if a is a cell, subdivide 
it into its Whitney decomposition of cubes. The sum converges in the 
mass norm to a chainlet J a . □ 
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Later, we will see that J\ uj — f a uj in the classical setting, where the 
LHS is the chainlet integral of JUand the RHS is the Riemann integral. 

A polyhedral chain is defined to be a chain of /c-cells. P = a « <T «- 
Every polyhedral chain is therefore a chainlet. 

THEOREM 5.14. Polyhedral chains are dense in chainlets. 

Proof. This follows since monopolar chains are dense in chainlets and 
polyhedral chains are limits of monopolar chains. □ 

Polyhedral chains include simplicial chains which are important in 
algebraic topology. Singular chains are also dense in chainlets. There 
is a natural chainlet Mayer- Vietoris result for chainlet homology using 
operators in this paper which will give a method for computing chainlet 
homology classes. 

5.5.3. Smooth submanifolds with boundary in R n . An oriented, smooth 
fc-submanifold N in M n may be triangulated. Each simplex <jj of the 
triangulation determines a fc-direction and a fc-mass. Let cti be the 
unique fc-cell with the same mass and fc-direction. Choose a point Pi 
in the simplex, or near it. The monopolar chain XXP«> a «) * s Cauchy 
in the 1-natural norm. The limit chainlet Jn represents N in that 
J n uj = J Jn uj for all forms uj G B\. (The LHS is the Riemann inte- 
gral over submanifolds, and the RHS is the chainlet integral.) Smooth 
submanifolds are dense in chainlets since polyhedral chains are dense. 

5.5.4. Differential forms. We embed B r k in j\f£. A /c-form uj is defined 
as a linear functional of monopolar chains. By the Riesz representation 
theorem, for each p G V, uj at p is represented by a /c-vector (3 such that 
uj(p; a) =< (3, a > for all simple fc-vectors a. Recall Pi = Yl 2~ nk {Pi] Q) 
converges to the improper chainlet V. Define Bi = 2~ nk {Pi] A) where 
A = P{pi)- Then Bi converges to an improper chainlet J w . Furthermore, 

Lemma 5.15. 



and \ Juj^ r = \u\^ r 
Proof. 

/ = lim / r] = lim 2~ nk r](pi; ft) = lim ^ 2~ nk r]A*uj(p i ) = I r]A*u. 

□ 



7] A-kUO 
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6. Measure theory and chainlets 

6.1. Lower semicontinuity. Define 

\u\l r = max{|H| ,...,p|M|}. 

Then the norms \uo\^ are decreasing as p — > to the norm ja;!^- 1 . It 
follows that the dual norms \P\^ are increasing to the norm |P|'' r_1 . 
The Banach spaces |P|p r are isomorphic to J\f£ for each p. (The spaces 
are the same, only the norms are different.) So each norm \P\ P r is con- 
tinuous in Ml- This proves that the limiting norm, the (r — l)-natural 
norm, is lower semi-continuous in _A/£. We similarly show that each 
\P\^ S is lower semicontinuous in Mjl for each < s < r, by multiplying 
the higher order terms by p. 
If J G JVfc and < s < r define 

|J| S = inf{liminf \P\^ : h J}. 

This quantity is the s-norm in the r-natural space. It coincides with 
the s-natural norm if J is a monopolar chain P. This follows from 
lower semicontinuity of | • | s . In particular, we may now freely speak of 
the mass of a chainlet J, realizing that the mass might be infinite. 

Theorem 6.1. If J e Nl then \T U J - J\^ < |«||J| r _i. 

Proof. Choose Pj — > J with [Pl^- 1 — > | J\ r _i. The result follows since 

\T U P-P\^ < |u||P| r _i. □ 

6.2. The operator Vec j . We define an operator Vec j : _A/£ — > A^.(V) 
as follows. It will be fundamental to integration, differentiation and 
measure of chainlet domains. 

For P = ^2(pi',o^), define 

Vej{P)=Y,°i- 

In what follows, set j — 0. The general case (proved in a sequel) is 
similar, but is not needed for this draft. 

Lemma 6.2. Suppose P is a monopolar chain and \\u)\\o < oo. Ifw(p) = 
u for a fixed coelement uiq and for all p, then 

J^co = uj (Vec (P)). 

Proof This follows for monopolar chains P G P° since u (a) = Uo(Vec°( 

□ 
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THEOREM 6.3. If P is a monopolar k-chain and r > 1, then 

M(Vec°(P)) <\P\^ r . 
If supp(P) C B £ (p) for some p G V and e > 0, then 
\P\h < M(Vec°(P))+eM(P). 

Proof Let rj be a coelement such that |?7o|o = 1, and rjo{Vec°(P)) = 
M(Vec°(P)). Define the fc-form rj by f]{<y p ) := r]o(a). Since rj is constant 
it follows that ||?7|| r = for all r > and ||<ir/|| r = for all r > 0. Hence 
\v\r — Mo — \Vo\o — 1- By Lemma f6. 21 and Theorem 15.61 it follows that 

M(Vec°(P)) = r] (Vec (P)) = f V < \v\r\P\^ = \ p \^- 
For the second inequality we use the definition of the r-natural norm. 

I f I 

It suffices to show that 4jw^ is less than or equal the right hand side 

for any 1-form uj of class B 1 . Given such uj define the /c-form uu (a q ) : = 
u)(a p ) for all q. By Lemma f6. 21 



Is 


< 




+ 


J 10 - UJ 













< \u(p)(Vec°(P))\ + sup \uj{p) - Lu(q)\M(P) 

qGsupp(P) 

< \\u\\ M(Vec°(P)) +e\\uj\\ 1 M(P) 

< \uj\^ (M(Vec°(P)) +sM(P)) 

□ 

Corollary 6.4. If J is a k-chainlet of class J\f r , r > 1, then 

M(Vec°(J)) < \ 

Proof. By Theorem 16.31 the result holds for monopolar chains. Choose 
Pi^ J with M(Pi) -> M(J). Then Vec°(Pi) -> Vec°(J) by continuity 
of Vec. By lower semicontinuity of mass, 

M(Vec°(J)) < liminf M(Vec°(Pi)) < hminf \P$ r = |J|K 

□ 

Theorem 6.5. Let Ji be a sequence of k-chainlets of class Af r such 
that for some C > and k-vector a the following hold: 

M(Ji) < C, supp(Ji) C B e .(p), Vec°(Ji) —> a as e € —> 0. 

Then J = lim Ji exists and Vec (J) = a. 
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Proof. Suppose Ej < S{ for j > i. Choose monopolar chains P s — > Ji—Jj 
with M(P S ) -> M(Ji - Jj) and supp{P s ) C B e .(p). By Theorem IO 

IP.I" 1 <M(Vec°(P s ))+eM(P s ). 

Hence 

\Ji - Jj\ kr < \ Ji - Jjf 1 < liminf M(Vec°(P s )) + - J,). 

But ^ec (P s ) — > Vec°(Jj — J^) — > by continuity of Kec. Hence 
J = lim J s exists and Vec°(J) = a. □ 

Corollary 6.6. Vec°(dJ) = 0. 

Proof. It suffices to prove this for (0, fc)-poles (3, by continuity of the 
operators. But Vec°(p; V u a) = implies Vec°{d(3) = 0. □ 

Given a monopolar chain P = Y^, a i a i, let AP denote its normalized 
linear contraction: AP = ^2(Xpf, \ k cti)/\ k . This contraction operator 
is continuous and extends to chainlets. By Theorem 16.51 Vec°( J) = 
lim A ^ A J. 

6.3. Part of a chainlet in a Borel set. The norm on a finitely 
additive fc-vector valued set function \i is defined as 



|/i| qr = sup |y u-dfi: |u;| qr < 1 

We can extend Whitney's theorem (jWj, XI Theorem 11 A) to as- 
sociate a fc-chainlet J with finite mass to a finitely additive fc-vector 
valued set function \ij defined by 

fij(X) = Vec°(J[ x ). 

Theorem 6.7 (Chainlet representation theorem). If J is a k-chainlet 
with finite mass, there exists a unique k-vector valued set function fij 
such that 



to = J u) ■ d/j. 

The correspondence is an isomorphism such that 

|/ij| q " = |J| q ". 

Whitney's proof extends, but the methods in this paper lead to a 
more direct proof. 

The ideas extend to define j-polar /c-vector valued set functions, 
leading to a similar representation theorem for chainlets with \J\j < oo. 

The isomorphism preserves operators. 
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6.4. Numerical method. Fix a mesh of V. Let J be a chainlet with 
finite mass in V. Let U,i be a cell in the mesh. We can define the 
part of J in Ui and denote it by J\vi- Then Vec°(J|t/J is a k- vector. 
Let Pi be a point in t/j. The sum X^ ec °(^VJ converges to J in the 
natural norm. That means we may do calculus over this monopolar 
chain and be assured that our results converge to the continuum limit. 
The operators involved in defining these monopolar chains all may be 
implemented numerically. 

6.5. Graphs of L 1 functions, (sketch) A nonnegative L 1 function 
/ : [0, 1] — > M is an increasing limit of simple functions gt each of which 
is a sum of indicator functions satisfying 



where /i denotes Lebesgue measure. The graph of an indicator function 
£s is a 1-chainlet with finite mass by Theorem 16.71 Denote the graph 
of a simple function g by T. The chainlet Tj — is a finite sum of 
difference chainlets. We estimate the total 1-norm of these difference 
chainlets. Since the Lebesgue area of the subgraph of / is finite, the 
1-norm of Tj — tends to zero as i, k — > oo. (Use Theorem 16.11 ) 



7.1. Multiplication by a function. A zero form e Bq is a function 
: V — > R. We may write 0(p) = (f)(p; 1), depending on the context. 
Recall |0|^ =max{||0|| o ,...,||0|| r }. 

Define J2(Pu «i) = J2{Pil 0(Pi)«i)- 
Proposition 7.1. (i) <f>(P + Q) = <j>{P) + 0(Q) 

(ii) (0 + i/))P = <pP + ^P 

(iii) {U)P = 4>&P) 

(iv) 0P = aP if 0(p) = a for all p 

(v) M((0 + <0)-P) = M(0P) + M(i)P) if (j)(p)^{p) > 

(vi) M(0P) < M{i)P) if0< <f)(p) < ip(p). 

Proof. These follow readily from the definitions. □ 

Proposition 7.2. (i) M(0<9P - <90P) < k\<j)\^M(P) 
(ii) M(d<pP) < k\<p\^M(P) + \<p\^M{dP). 

Proof. These are essentially the same as in jW], p. 209. □ 

Theorem 7.3. // e B r Q and P e V k , then 




7. Vector fields and chainlets 




r 
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Proof. Let e > 0. By Theorem 12 A\ there exists P = 5^j = o ^ sucn that 
G 7^ and |P|^ > Yl \\ Dj \\j - £ - lt suffices to prove that 



i=0 

Clearly, this holds for j = 0. Assume it holds for norms less than j. 
Let U = (u, U') and = A^V. Then 

A{ja = A^A^a) = A u (3. 

Therefore, 

\<j>(p; A^ = \<j>(p; A u f3)\^ 

< \<j)T u (p-(3) -T u 0(p;/3)|^ + |T u 0(p;/3) -0(p;/3)|^. 

Now 

\<f>T u {p^)-T^p)\^ = mp + u)-4>{j>)){p + u;(3)\^ 

^EpTMi^-^i^ 

i=0 ^ ' 
3-1 



i=0 ^ ' 

^ E ( ' 7 1 ) l*M|A{ra|| 



-1 



i=0 

and 

it u #(-#^ < iriii^i^- 1 < E ( J ; 1 ) i^mi/^- 1 

i-i 



=o 

It follows that 

j 



n V / 



|0(p;A^)|^<^(^)|0|^||A^|| r 

i=0 ^ ' 

Define (f>u){p; a) = u(p; <p(p)a). It follows by linearity that 

(j)Uj(P) = u(<pP). 
The next result follows from Proposition 17.11 



□ 
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Proposition 7.4. (i) <p(uj + rj) = <j){u) + 0(77) 

(ii) (4> + ip)uj = (pu + ipoj 

(iii) = 0(^) 

(iv) 0u; = au; z/ (j){p) = a for all p 

(v) |(0 + <0Vl te = 1 1 1,0 + if(j)(p),il>(p)>0 

(vi) < |^cj| fe z/0 < 0(p) < V»(p) 

The next two propositions follow from Proposition 17.21 and Theorem 

Proposition 7.5. (i) \4>duj - dcpu^ < (k + l)|0| 1,1 |a;| te 
(ii) \d(puj\^ <{k + + \<p\^\du\*°. 

Proposition 7.6. 

l<M* r <E( • Wi^- 

7.2. ^-vector fields. A ^-vector field X on V is defined to be a func- 
tion X : V — ■> Afe(V). Therefore, functions are a field of monopolar 
0-vectors. We next consider monopolar fc-vector fields. 

The vector field e± in V is associated to the 1-chainlet J ei defined 
as follows: We first define J ei in a unit ra-cube Q. Let Q = R x I 
where R is an (n — l)-face of Q and / is a 1-cell with direction e\. 
Let po be the midpoint of R. Subdivide R into binary (n — l)-cubes 
Rk,i with side length 2~ fc . Denote the midpoints of the Rk,i by pk,i- Let 
P = °"o denote the 1-cell supported in Q, with direction e±, and with 
one endpoint p . Let Pk = ^2 ^~ nk Ck,i with each 07^ parallel to Oq and 
with endpoint p^. Then M(Pjt) = 1 and the sequence {P&} is Cauchy 
in the 1-norm since the distance of translation tends to zero and the 
total mass that is translated i M(Q). (See Theorem 12.41 or |H7j for 
details.) Denote the limit by J ei . Use the preferred basis to construct 
n-chainlets J ei associated to the vector field e*. 

A vector field X can be written X = ^™ =1 0ie«. Therefore, the 
chainlet Jx = Y^=i * s associated to X. Define 

n 

\x\^ = Hit- 

i=l 

7.3. Exterior product ExJ- Let (3 G A m (V). Define the exterior 
product operator 

■ 1 k ' k+rn 

by 

The next two results follow directly from the definitions. 
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Lemma 7.7. EpP is a bilinear operator in v and P satifsying 

(i) EpEp = 

(ii) Ep(P -Q)=EpP-Q + (-l) k P ■ EpQ 

(iii) UEpP = E f .pf*P 

(iv) <f>EpP = Ep<j)P = E^pP 

(v) \EpP\^ = |Ppr. 

Proposition 7.8. 

< M((3)\P\^. 

If / : M — > N is a mapping and Y is /-related to X, that is, /* oX = 
Y o /, then Ey(f*P) = f*(ExP)- Of course, is / is a diffeomorphism, 
then /*X is always /-related to X. 

Define the interior product operator of forms by duality 

as 

ipuj(p; a) = u(p;Epa). 

Lemma 17.71 yields 

Lemma 7.9. ipu is a bilinear operator in v and u satisfying 

(i) ipip = 

(ii) ip{uj A 77) = ipuj A 77 + (— l) k u> A iprj 

(iii) /*(^w) =i ft pf*uj 

(iv) = ^^u; = i^cj; I 

(v) |2/3Co>| 1,r = |cup r . 

THEOREM 7.10. Let X be an m-vector field on V, and P a monopolar 
k-chain. There exist constants C(X,k,m,r) > such that 

\E x {P)t < C(X,k,m,r)\X\^\P\^. 

Proof. Set X = J2<f>i e i- Observe that \E ei (P)\^ < \P\^. Then 

"(E^P) = ieMP) = ieMP) < sup ^CP) = |p| hr _ 

Therefore, by Lemma 17.71 (iv) and Theorem 17.31 

\E^(P)\^ = \<P(E ei P)\^ < CtfarM^Pfr 

where the constants Cj(0, r) = J2i=o y \ J l&l^- Now take sums over 
the basis (ex, ... , e n ). To obtain the result, set C(0, r) = ^ Ci{<P, r )- 

□ 
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If J is a chainlet in Afk r , X and m- vector field, and Pi — > J , define 
E x J = ]imE x P i eAf£ +m . 
Lemma 17.71 (iii) extends to J: 

f*E x J = Ef t xf*J 

and 

f*i x uj = if,xf*u- 
A similar construction leads to 

E X iJ 

for a j-polar m- vector field XK As with most of the operators, j = 
1,2, ... , we obtain fields polypolar vectors expressed as chainlets. 

Since the operators and products commute with pushforward, the 
definitions and relations extend to smooth manifolds. 

If / : M — > N is a smooth mapping and Y is /-related to X, then 
Y = f*X and i Y f*uj = f*i x co. 

7.4. Lie derivative C X J- If X is a vector field and a is a k- vector 
field in a manifold, we recall the classical Lie derivative of a in the 
direction X. Let g t be the flow of X. Fix xq and define 



£x(ol) xo := lim t . 



t 

If X is smooth, then the limit exists. 

We similarly define C x on fields of (j, fc)-poles The limit exists in 
the j-natural norm as a field of (j, fc)-poles. This leads to a definition 
of C X J for a chainlet J. 

Lemma 7.11. 

C x d = dC X - 



Proof. 



f- t da Xt - da X0 f- t a Xt - a xo 



C x da = lim J ~ L xt x ° = d lim J ~ l xt x ° = dC x a. 
t t 

□ 

Define 

£ x u(p; a) := u(p; C x a). 
C x ui is bilinear in X, uj and C x uj A T] = C x uj A rj + u> A C x r]. 

Theorem 7.12. (i) Vec°{V x ) = C x 
(ii) dV x = V x d. 
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Proof. The first part follows by continuity of the operator Vec. 
d£xu(p; a) = Cxu(p] da) = tj(p; Cxda) = u(p; dCxCt) = Cxdw(p; a). 
The second follows from the definitions. □ 
Theorem 7.13. V x = dE x + E x d. 

Corollary 7.14 (Cartan's magic formula). C x u = di x oj + ixduj. 

Proof. By Theorem ITT21 £vP = Vec°(V x P) = Vec°(dE x P + E x dP). 
It follows that dixUj(P) + i x du(P) = uo(d(P A X)) + uj(8P A X) = 
u(V x P) = Cxuj(P). □ 

For a diffeomorphism /, 

f*C x uj = C JmX f*u. 

If / : M — * N is a mapping and Y is /-related to X, then 

£ Y f*uj = f*C x u. 

7.5. Translation T X J . The translation operator T u (p; a) = (p + u; a) 
is continuous in the chainlet norm. Similar techniques as those given 
above show that this operator extends to TxJ for smooth vector fields 
X and chainlets J. 

8. The chainlet complex in a manifold 

Now set V = E n . A singular k-chainlet gj in a smooth n-manifold 
M of class M r is a continuous mapping g : R n — ► M and a /c-chainlet 
J in R n . Consider the vector space J\f£(M), generated by all singular 
fc-chainlets subject to the relation: g(p; a) = h(q; (3) if and only (h^ 1 o 
g)*(g(p; a)) = h(q; /3), for all (p; a) £ 1" x A°. The boundary operator 
is well defined on chainlets in a manifold by Corollary 13.31 We obtain 
a chain complex of vector spaces since d o d = 0. 

K +r (M) 4 • • • *> M£}(M) ± N r h {M) ± NEliM) ± ■ ■ ■ ^ A/T r (M). 

Forms of class B r k and of class CI are well defined on M. The relation 
f*u = lo o ft, leads to a well-defined integral of smooth forms of class 
B r k over chainlets in M k {M). 




The operators V„ and d are well-defined on manifolds since they 
commute with pushforward of diffeomorphisms. 

We deduce extensions of Theorem 15.71 and Stokes' theorem 15.81 
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Theorem 8.1. 



D u uj 



10 



for u G B r k (M) and J G M r k +l . 

Theorem 8.2 (Stokes' theorem for chainlets in a smooth manifold). 

duj = 



10 



for lu G B r k (M) and J G Af^. 



In Riemannian manifolds, the translation operator T u is defined via 
the covariant derivative. The operators _L and Hodge star are well 
defined, leading to a divergence theorem for chainlets in Riemannian 
manifolds. 

9. Integration of rough forms over rough domains 
Given a chainlet J G Ml and a form uo G Bf., < s < r, define 

h J}. 



LO 



Then 



to 



du 



inf 






j Pi 




< 


h 




JdJ 



J d-k uo 


< 


I - 






J±dJ 



and 



J D u u 


< 


I ^ 






JVuJ 



< \u\\J\.\oj\ 



Whitney's example |W2j of a function nonconstant on a connected 
set of critical points shows this inequality is sharp. (See also J\] . ) 

10. Further operators 

10.1. Slant product J/X. Suppose a G A° k and (3 ^ G with 
< m < k < n. Define the slant product 

Sp a = a/0 := (-l) k ^ ± ((3A ± a)/M(/3) 2 G A° k _ m . 

If m = 0, slant product reduces to division of a fc-vector by a nonzero 
scalar. Therefore, if k = m = 0, slant product reduces to division of 
real numbers in R 1 . 
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Lemma 10.1. The slant product Sp : A° — > A°_ m is linear and satisfies 

(i) Vec°((3) c Vec°(a) (3 A (a//3) = a 

(ii) \/ec (/3) JL Vec°{a) =^ {a /\ (3) / (3 = a and a/ (3 = G F 

(iii) a/a = 1. 

Proposition 10.2. \P/(3\^ < \P\^/M((3). 
Proof. Then 

\P//3\* r = |/3A 1 P|^/M(/5) 2 < |P| h 7M(/3). 

□ 

Therefore, the slant operator extends to chainlets in E n . 
\J/P\\ kr < \J\^/M((3). 

Lemma 10.3. Slant product is a linear transformation 

Sp-'Pl^ ^k-rn 

satisfying 

(i) SpSp = 

(ii) f*oS p = 5/,j8 o /* 

(iii) (frSp = Spcj). 

Therefore, slant product extends to chainlets in manifolds. 

If X is a nonzero monopolar fc-vector field, slant product is defined 
at each point. Let us further assume that the mass of each k- vectors 
is bounded below by a constant K. 

Theorem 10.4. If X G B r , there exists a constant C(r,X) such that 

\P/X\^ < C(r,X)\X\^\P\^. 
Proof. By Theorem 17.31 

| JL ((0ei)A J_ P)\^ = |0(eiA JL P)\^ < C(r, 0)|P|K 
Suppose X = *YTi=\$i e i- Let C(r,X) = ^C(r, Xj). Then 
\P/Xfr < I J- (<PMPs)A J- Pt/M{X{p s )f < K 2 C(r,X)\X\^\P\^ 

□ 

It follows that J/X is well defined for a nonzero k- vector field X. 
Order j vector fields are treated in a similar fashion. 

For each (3 G A^, define the extrusion operator on forms 

extp : B r k -> B r k+m 

extpui{p; a) = u(p; a/ (3). 
This extends to fc-vector fields X. 
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ext x uj(p; a) = u(p; a/X(p)). 

Define Tip = Sextp + extpS where 5 = *d * . It is an open question 
whether or not Tip is a derivation. 

Slant product is not generally defined on smooth manifolds since 
pushforward does not commute with the perp operator. 

Cross product J x X. Define 
by 

a x (3 :—± (a A (3). 

For ki = k 2 = n/3 this product combines pairs of (0, fc)-poles and 
produces a (0, fc)-pole. Of course, when n = 3 this corresponds to the 
standard cross product of vectors. 

In the degenerate case with v = w, then Vec°(v) A Vec°(w) is zero. 
Thus _L of it is zero. On the other hand, if v is orthogonal to w, then 
the wedge product is a (0, 2)-pole with mass the same as the product 
|i>||w|, so its perp corresponds to a vector with norm \v\\w\ that is 
orthogonal to this (0, 2)-pole. 

We next show this extends to J x X where X is a smooth fc-vector 
field in V. 

Theorem 10.5. 

\P x X\^ < C(r,X)\P\^\X\^. 
Proof. Let X = tfiiCi. Then by definition and Theorem 17.31 

\P x (j)^ = ^ | _L ((p s ,a 8 ) A (j)i(p s )ei)\^ 
= ^\HPs){Ps,a s A ei)\^ 
<C(r, 0010^1^1^- 
The result follows by setting C(r, X) = C(r, (f>i). 

□ 



As usual, we define a dual operator on forms as 
(ijj x f3)(p; a) = u(p; a x j3). 
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Intersection product JClX. Let ki + k 2 > n. Let a G Vl\, j3 G V 3 ? . 
Define 

an (3 :=_L (_L aA _L (3). 
For ji = j2 = this identifies the intersection product of a and 

As above, we may extend this to J H X where X is a fc-vector field. 
The dual operator on forms is define by 

uj fl X(p; a) = uj{p; a fl X). 

Projection. If a G Vk+j, (3 G Vk define 

TTpot :=_L (_L a n /?) n /?. 

As above, we may extend this to nxJ where X is a k- vector field. The 
dual operator on forms is define by 

ti x uj(p; a) = uj(p; TT X a). 

11. Applications 

11.1. Solutions to Plateau's Problem. Suppose a is a 2-cell in 3- 

space. Let u be a unit vector orthogonal to Vec°(cr). Then V u a is a 
dipole surface that locally models a soap film without branches. By 
adding three of these along a common dipole edge at angles of 120 dog , 
we can obtain a branched surface. 

A mass cell is defined to be e u a. Then de u a = V u cr + e u (da). Sums 
of mass cells and dipole cells give models of soap films with curvature. 

We call a chain of dipole fc-cells and mass fc-cells a fc-dipolyhedron. 

As an application of chainlet methods, we observe that any naturally 
arising soap film S spanning a smooth curve 7 can be expressed as a 
limit of dipolyhedra S = lim Di in the natural norm with ODi supported 
in 7. The boundary of S is also supported in 7. 

A similar construction in R 4 |H8j leads to existence of minimal span- 
ning set for a fixed Jordan curve, a solution to Plateau's problem, 
assuming a bound on energy. The minimizer has soap film regularity 
a.e. jH.9. and has surface area smaller than any other Plateau solutions 
to date. 

The new methods of this paper shed some light on the general prob- 
lem which is under investigation. 

11.2. Chainlets and distributions. The question arises whether or 
not we may replace currents with the smaller space of chainlets in 
analysis. The simplest case is k = 0. In this section we ask whether we 
can "take the derivative" of 0-chainlets. Distributions are defined over 
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test functions with compact support. For O-chainlets, test functions 
need only be integrable, in order to take the derivative 

If J is a O-chainlet in R, define J' = V u J where u is the unit vector e\. 
This corresponds to the derivative of a smooth function /. A smooth 
function / determines a chainlet Jf if its integral is finite. 

Lemma 11.1. V U J/ = -Jf. 

Proof. The proof reduces to showing 

/ » = -/» 

for all smooth g. That is, 




We use integration by parts. f(b)g(b) — f(a)g(a) = f fg' + J f'g. Since 
the integral of / is finite and / is smooth, the left hand side tends to 
zero for a — > oo and b — > — oo. □ 

Chainlets offer more structure, and therefore a richer theory that is 
easier to work with than distributions since they form a normed space. 
Applications to PDE's are anticipated. 

11.3. Dual mesh convergence. In the discrete theory of simplicial 

complexes, the problem of existence and convergence of a geometric 
Hodge star operator has eluded mathematicians until recently. Barycen- 
tric subdivisions or or barycentric duals are commonly used. The cir- 
cumcentric dual is introduced in |D, H, L, M| which has the advantage 
of orthogonality to the primary mesh. But convergence to the smooth 
continuum remained elusive. The author proved the chainlet discrete 
Hodge star converged to the smooth continuum QH6J, see also |H7J). 
Wilson WiJ constructed a combinatorial Hodge star operator via a 
dual mesh, and showed convergence to the continuum as the mesh of 
a triangulation tends to zero. The methods of |H6j apply to simpli- 
cial complexes as a special case, but there is no mention there of dual 
complexes. We address this now. 

Let (S, T) be an ordered pair of two simplicial complexes. We will 
call (S, T) a mesh pair if the fc-simplexes o\. of S are in 1 — 1 corre- 
spondence with the (n— fc)-simplexes r n _^ of T. We call S the primary 
mesh and T the dual mesh. Let pk be the barycenter of and q n -k 
the barycenter of r n _jfc. Let a n ^k =-L Vec°(o~k) and j3 n -k = Vec°(r n ^k). 
We say that a sequence of mesh pairs (S l , T l ) is a Hodge sequence if 

(i) The mesh sizes of both S % and T % tend to zero as i — ► oo. 



CHAINLETS 



45 



M(fa_ k - n ■ 
-—7-r- — : > U, as i — > oo 

(iii) 

\v\~<in-k\ -^°> asi^oo. 
LEMMA 11.2. If (S\T l ) is a Hodge sequence then 
It* - 1 at 



as i — > oo. 



MK_ k ) 
Proof. By the triangle inequality, 

+ \Tp k - qn _ k T n - k - (pfcj/Jn-fc)!" 1 

+ |(Pk,A»-fc) - (Pk,a n - k )\ h 
+ \(Pk,ain-k)- -L crfcl 111 



We know 



\(ph<-k)- - 1 4 



-> 



by the definition of the operator J_ . In particular, the distance of 
translation of estimates of (p\, at l n _ k ) — J_ o\ is tending to zero, the 
total mass is kept constant. (See Theorem 12.41 ) The second term is 
similar. 
Also, 

\ T% n-k ~ r ^P i k -q i n _ k T n-k\^ — \Pk ~ Qn-k\-^( T n-k)- 

The third term simplifies to 

□ 



Since _L is a continuous chainlet operator and dual to Hodge star 
* on forms, this lemma guarantees that a Hodge sequence will limit 
to the smooth continuum w.r.t. the Hodge star and boundary opera- 
tors. Therefore, discrete methods relying on such meshes will produce 
reliable discrete approximations to Gauss and Stokes' theorems. Nu- 
merical estimates on convergence may be gleaned from the definition 
and lemma. 
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11.4. Wavelets over chainlets. A chainlet transform is defined 
on R n . (Sketch) Cover R n with a binary grid Q. Let J be a fc-chainlet. 
For each Q m in the grid, the part of J in Q m is a chainlet J[_Q m - 
Since mass is lower semi-continuous in the chainlet norm |H6j . a m = 
M(J[g m ) is well-defined. Let p m denote the midpoint of Q m . Let 
T( J, Q) = Yl a m Vec°(J[Q m ). Then T( J, Q) is a sequence of monopolar 
chains converging to J in the 1-natural norm. 

By setting P = T(J,Q ) and P { = T(J,Q i+1 ) - T(J,Qi), we may 
write 

oo 
i=l 

The Pj may be written in terms of an orthonormal basis of the sub- 
space of monopolar chains Pj = Ys=i so that 

OO fc!(n — fc)! 

i=l s=l 

Wavelets over chainlets now become natural, as do Fourier series over 
chainlets. If u = Yl'jLi Vj is a wavelet series or Fourier series, then 

oo oo „ 

^■ = EE / 

i=l j = l JP i 

We end up with a convergent double series of trivial integrals. All of the 
operators involved may be implemented numerically since monopolar 
chains are dense in chainlets. 




Appendix A 

11.5. Deduction of the coordinate calculus. We next see that 
standard integral equations of coordinate calculus follow readily from 
the concise integral equations (i)-(iii) of chainlets. This paper, there- 
fore, provides the foundations for a full theory of calculus starting from 
basic principles of multilinear algebra. § The following proofs are not 
possible from the perspective of the Cartan theory which is missing the 
prederivative operator V u and the geometric Hodge star operator _L. 
We first retrieve the standard formulation of the exterior derivative of 
a multivariable function / : E n — >• R: Let (ei, . . . , e n ) denote the stan- 
dard basis of R n . The pair (x; v) denotes a 1-vector v in the tangent 

§ Indeed, in the spring of 2006, the author will be teaching an experimental course 
on chainlet theory to 30 math students at Berkeley, all of whom have studied 
abstract linear algebra, with no assumptions of calculus, algebra, or real analysis. 
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space of i G R". A partial derivative of / at x G V is defined as a 
directional derivative 

^-(x) = D et f(x;l) 
where (ei, . . . , e n ) is a preferred basis and x = x\e\ + • • • x n e n . 

£-< 

oxi 



Theorem 11.3. df(x;ei) = ^-dxi(x;ei 



Proof. 

df(x; = f(x; de*) = f(x; V ei (l)) = D e J(x; 1). 

Corollary 11.4. df = ^ §^^ x i- 

The Fundamental Theorem of Calculus follows from (i): 

f f'{x)dx = [ df= [ f = I f = f(b) - f(a). 

J a J\aM Jd\aM J \b\-M 



□ 



Green's theorem over a bounded open set (/CI also follows from (i): 

f f f dQ dP 

/ Pdx + Qdy = / d(Pdx + Qdy) = / ——dxdy — —dxdy. 
Jdu Ju ' Ju ox dy 

The chainlet divergence theorem follows from (i) and (iii): 

/ d-k UJ = I UJ. 
J J J±dJ 

For a simple domain V in R 3 , we retrieve the Gauss' divergence theorem 
for 1-forms F defined in a neighborhood of V. Suppose F = Pdx + 
Qdy + Rdz. Then 

d -k F = d-k (Pdx + Qdy + Rdz) = d(Pdydz + Qdzdx + Rdxdy) 
OP dQ dR\ , , , „ ^ JT , 
ax ay az / 



Thus 



V-FdV= d*F= F= F-ndS. 

V JV J±dV JdV 



The chainlet curl theorem takes the form Jj-kduj = J 9±j uj. Therefore, 
if S is a surface in M 3 and F a smooth 1-form defined in a neighborhood 
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of S in K 3 we calculate 

*dF = *d(Pdx + Qdy + Rdz) 

fdR _dQ\ dx+ fdP _ dR\ ^ + fdQ _ dP\ ^ 
\dy dz J \dz dx J \dx dy J 
= V x F 

Therefore, 



/ (V x F) ■ ndA = [ *dF= [ F = I F ■ ds. 

Js J±S JdS JdS 
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Appendix B 

11.6. Relation to previous versions. Chainlets have been develop- 
ing over a period of years. In earlier versions, they relied on coordinate 
expressions which made proofs long and relatively cumbersome. They 
initially took the viewpoint of Whitney, starting with completions of 
polyhedral chains with respect to a norm. The Ravello lecture notes 
|H6j started with polyhedral chains and developed the first discrete 
chainlet theory with the definition of a k-element as a geometric rep- 
resentation of a k- vector, based at a point, now called a monopolar 
k-vector. A monopolar k-chain is a finite sum of monopolar /c-vectors 
supported in finitely many points, and these form a dense subset of 
chainlets, as seen in |H6j . The main object of the Berkeley lecture 
notes |H7j was to produce a coordinate free theory of chainlet calculus 
without any assumption of the integral theorems of classical calculus. 

The initial goal of this paper was to develop the theory without 
using polyhedral chains, except as examples. The Koszul complex aided 
enormously in this quest. Setting aside polyhedral chains soon led to 
large simplifications of the theory. Now there is essentially one classical 
limit to establish for the full calculus, beyond the Cauchy sequences 
needed for examples, namely that of Lemma 12.51 proving the existence 
of the prederivative operator. 



Conclusion 

The prederivative operator V u is an operator at the very foundations 
of mathematics. For V n leads to boundary, which, in turn leads to 
exterior derivative through duality. It links boundary d to exterior 
product e u (a) = (a A u) in a "magic formula" 

VxJ = de x J + e x dJ. 

The operator V M leads to the existence of norms on chainlets which, in 
turn, leads to uniform convergence theorems, not available for distribu- 
tions and currents. The prederivative, combined with the preintegral, 
leads to a general divergence theorem with three integrals, rather than 
two, that models well conservation of matter and energy in curved 
space. The full continuum with its polypoles of all dimensions and 
orders reveals mysteries of mathematics and physics. Both simple and 
multilayered complex systems naturally arise from this which are self- 
organizing. The polypoles are natural candidates for the basic compo- 
nents of a "genetic code" of the universe itself. 
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